Local field theory for disordered itinerant quantum ferromagnets 
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An effective field theory is derived that describes the quantum critical behavior of itinerant ferro- 
magnets in the presence of quenched disorder. In contrast to previous approaches, all soft modes 
are kept explicitly. The resulting effective theory is local and allows for an explicit perturbative 
treatment. It is shown that previous suggestions for the critical fixed point and the critical behavior 
are recovered under certain assumptions. The validity of these assumptions is discussed in the light 
of the existence of two different time scales. It is shown that, in contrast to previous suggestions, 
the correct fixed point action is not Gaussian, and that the previously proposed critical behavior 
was correct only up to logarithmic corrections. The connection with other theories of disordered 
interacting electrons, and in particular with the resolution of the runaway flow problem encountered 
in these theories, is also discussed. 
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I. INTRODUCTION 

The theory of quantum phase transitions, i.e. phase 
transitions at zero temperature (T = 0), is an important 
problem that is, relevant to many topics in condensed 
matter physics.El Perhaps the most obvious example is 
the T = transition from a paramagnet to an itinerant 
ferromagnet as it occurs in, e.g., diluted Ni, or in soHd 
solutions like MnSi. Historically, this was the first quaru 
turn phase transition that was studied in detail. Hertzcl 
showed how to treat this transition by means of renormal- 
ization group (RG) methods, and he concluded that the 
transition in the physically interesting dimension d — Z 
was mean-field like. This conclusion hinged on the obser- 
vation that in quantum statistical mechanics the effective 
dimension of a system for scaling purposes is d-\- z, with 
d the spatial dimensionality and z the dynamical criti- 
cal exponent. Since in a simple theory at tree level one 
has z = 3 for clean itinerant, quantum ferromagnets, and 
z = 4 for disordered onesj3 this seemed to imply that 
the upper critical dimension o?jf, above which one finds 
mean- field critical behavior, is d^ — 1 and (i+ = in the 
clean and disordered cases, respectively. 



This conclusion was later challenged. ETEI Reference ^ 
noted that Hertz's results for clean systems in rf = 1 — e 
dimensions were inconsistent with general scaling argu- 
ments. This left open the possibility of mean-field be- 
havior in physical dimensions. However, in Refs. ^js] it 
was shown that the critical behavior in d > 1 (clean case) 
and d > (disordered case), respectively, is not mean- 
field like after all. The salient point is that in itinerant 
electron systems at T = 0, soft modes other than the 
order parameter fluctuations exist. These modes are dif- 
fusive (in disordered systems) or ballistic (in clean ones) 
particle-hole excitations. Since they couple to the order 
parameter fluctuations, they influence the critical behav- 
ior. Specifically, they lead to an effective long-range in- 
teraction between the order parameter fluctuations. As 
a result, Refs. ^js] found that the critical behavior is gov- 
erned by a Gaussian flxed point which, however, does not 
yield mean-field exponents. The critical behavior deter- 
mined in these references was claimed to be exact. 

A separate, and seemingly unconnected, develop- 
ment in the many-electron problem has been the study 
of metal-insulator transitions of disordered interacting 
electrons. □ For one of the universality classes that oc- 
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cur in this problem, a transition was found that is 
not a r»etal-insulator transition, but rather of magnetic 
nature.l3 While the order parameter, and the nature of the 
ordered phase, could not be identified with the methods 
employed, the critical behavior for all quantities other 
than the order parameter was determined^ Apart from 
differences in logarithmic corrections to power laws, this 
critical behavior turned out to be identical with the 
Gaussian critical behavior for the disordered ferromag- 
netic transition. This led, in Ref. U, to the suggestion 
that the unidentified transition studied in Ref. |^ was the 
ferromagnetic transition. The discrepancy with respect 
to the logarithmic terms was explained as due to the fact 
that of the two integral equations in Ref. ^ only one had 
been shown to be exact. The proposal thus was that the 
two approaches describe the same transition, and that 
the critical behavior found in Ref. || was exact while the 
one in Ref. ^ represented an approximation. 

The theory developed in Refs. suffers from one ma- 
jor drawback: Since the additional soft modes were in- 
tegrated out in order to obtain a description entirely in 
terms of the order parameter, the effective field theory 
that was derived is nonlocafl and not very suitable for 
perturbatively calculating effects that depend on all of 
the soft modes in the system. The analysis in Refs. 
therefore was restricted to power counting arguments at 
tree level to show that all non-Gaussian terms are irrel- 
evant in a RG sense. While this turned out to be true, 
relying entirely on tree-level power counting can be dan- 
gerous. Indeed, even a one-loop analysis of Hertz's action 
would have revealed that the mean-field fixed point is un- 
stable, and the absence of such a calculation led to the 
instability to not be noticed for 20 years.t£l Furthermore, 
integrating out the fermionic degrees of freedom obscures 
the fact that the problem contains two time scales, a dif- 
fusive one and a critical one, which in itself makes power 
counting very subtle. This, combined with the sugges- 
tive relation between the ferromagnetic transition and 
the unidentified transition discussed above, and the puz- 
zling logarithmic discrepancies between the critical be- 
haviors found for the two transitions, makes it desirable 
to have a theoretical description of the quantum ferro- 
magnetic transition that takes the form of a local field 
theory which facilitates a controlled loop expansion and 
keeps the two time scales explicitly. Another motivation 
for constructing a local field theory is that it will allow for 
a more explicit study of the effects of rare regionsllil than 
was possible within the framework of Refs. |j,|5|, although 
we will not pursue this issue in the present paper. 

It is the purpose of thjC present paper, and a second one 
to be referred to as II,t3 to put these remaining questions 
to rest. We will focus on the disordered case, although 
we expect analogous conclusions to hold for clean sys- 
tems. By using a local field theory description, we will 
show that Ref. ^ missed effects of the two time scales 
that lead to logarithmic corrections to the Gaussian crit- 
ical behavior. Moreover, taking these effects into account 
leads to integral equations for the relevant vertex func- 
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tions that are identical to the ones derived in Ref. |. The 
current formulation makes it obvious that the transition 
described by these equations is the quantum ferromag- 
netic one, and it allows to determine the exponents in 
the ferromagnetic phase as well as those in the paramag- 
netic one. It furthermore shows that the integral equa- 
tions that were first derived in Ref. ^ are exact, and it 
elucidates many physical points that were rather obscure 
in Ref. ^ and to an even larger extent in Ref. 0. 

This paper is organized as follows. In Sec. jEJ we use 
methods developed in Ref. |l^ to derive an effective the- 
ory for disordered itinerant quantum ferromagncts that 
systematically separates massive modes from soft ones, 
and explicitly keeps all of the latter. In Sec. II] we give 
a RG analysis of this model. We first show how Hertz's 
fixed point, as well as the Gaussian fixed point of Ref. 

emerge within this framework. We then show that 
Hertz's fixed point is unstable, and the Gaussian one 
marginally unstable, due to the existence of two sepa- 
rate time scales, viz. a critical time scale associated with 
the order parameter fluctuations, and a diffusive one as- 
sociated with the additional particle- hole excitations. We 
identify an effective action that contains a stable critical 
fixed point. This action is not Gaussian, and its solu- 
tion is therefore nontrivial and deferred to II. In Sec. IV 
we discuss our results, in particular the relation of the 
present approach to previous ones, and the complications 
that the presence of two time scales leads to in scaling 
considerations. 



II. EFFECTIVE FIELD THEORY 

In this section we start with a simple model for inter- 
acting electrons in a disordered environment. We then 
introduce the ferromagnetic order parameter and iden- 
tify all other soft modes. Integrating out the massive 
modes leads to an effective field theory that describes all 
of the soft modes in the system. The general method 
employed here is the one that was developed in Ref. O. 



A. Model of itinerant electrons 

Our starting point is a general field theoretic represen- 
tation of the partition function of a ruajuy-fermion sys- 
tem, which can be written in the form,L3 



Z= / D[^,,p] exp(5[^,V]) . 



(2.1a) 



Here the functional integration measure is defined with 
respect to Grassmannian, or anticommuting, fields ip and 
'0, and S is the action. 



^ = ~j dr J dx ^aiy^,T) -^^pa{y^,T) - J dr H (t) . 

(2.1b) 
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We denote the spatial position by x, and the imagi- 
nary time by r. H{t) is the Hamiltonian in imaginary 
time representation, jS — l/T is the inverse temperature, 
a = 1, 2 denotes spin labels, and a summation over re- 
peated spin indices is implied. We choose units such that 
kB ^ h = = 1. The Hamiltonian describes a fluid of 
interacting electrons moving in a static random potential 



H{t) = J dy 



1 

2m 



VV'a(x, r) • V?/'a(x,r) 



-t-[u(x) - ^] -00 (x, t) V'a(x, r) 



x0a(x,r) . (2.2a) 

Here m is the electron mass, fi is the chemical potential, 
and m(x — y) is the electron-electron interaction poten- 
tial. For simplicity, we assume that the random potential 
w(x) is delta-correlated and obeys a Gaussian distribu- 
tion P[w(x)] with second moment 



{u(x)w(y)}ji, 



1 



2ttNfTc 



^(x-y) 



where 



{...}dis = J D[v]P[v]i...) , 



(2.2b) 



(2.2c) 



denotes the disorder average, Np is the bare density of 
states per spin at the Fermi level, and t^i is the bare 
electron elastic mean-free time. Our results will not be 
sensitive to the simplifications inherent in the assump- 
tions that lead to Eq. ( 2.2b ). We also mention that it 
would be possible to include more realistic features, e.g. 
band structure, in the model. However, ultimately we 
will be interested in universal behavior at a phase tran- 
sition that is independent of all microscopic details. For 
our purposes it thcrciojrc is sufficient to study the model 
defined in Eqs. (p|)£3 

As in Ref. |, and as is standard practice in the theory 
of magnetism, we break the interaction part of the action 
5, which we denote by S'int, into spin-singlet and spin- 
triplet contributions, S* For simplicity, we assume 
that the interactions are short-ranged in both of these 



channels .El The spin-triplet interaction S-^^l describes in- 
teractions between spin-density fluctuations. This is the 
interaction that causes ferromagnetism, and it therefore 
needs to be considered separately. We thus write 



(t) 



S = Sn + S, 



it) 



with 



c(*) _ 
•^int — 



Ft 
2 



J dx dr Us (x, t) • (x, - 



(2.3) 



(2.4a) 



where ris is the electron spin-density vector with compo- 
nents, 



'^s (x, t) = (x, r) erf 06 (x, t) 



(2.4b) 



Here the ai (i = 1, 2, 3) are the Pauli matrices, and Ft is 
the spin-triplet interaction amplit ude t hat is related to 
the interaction potential u in Eq. ( 2.2a ) via 



Ff = ^ / dx u(x) 



(2.4c) 



5*0 in Eq. (^^) contains all other contributions to the 
action. It reads explicitly. 



So = 



/3 



dr / rfx 



d 

V'a(x,T) — ^a(x,T) 



-■0a(x,r) \^—+flj-ipa{x,T) 

'-^ J dr J dx nc{yi,T)nc{x,T) 



dr J dx v{x.)lpai'^,T)lpa{'^,T) , 

with Tie the electron charge or number density, 
"-c(x, r) = ?/!q(x, r)?AQ(x, r) , 



(2.5a) 



(2.5b) 



and Ts the spin-singlet interaction amplitude. 

Before we proceed, we integrate put the quenched dis- 
order by means of the replica trick. EZI Pe rforming the dis- 
order average as prescribed in Eq. ( 2.2c ) replaces the last 
contribution to the action 5o, Eq. ( ^.5a| ), by 



N 

E 



47rA^Froi — jq 

ai,a2 — 1 



drdr' / dx C(x,-r) 



x0r(x,r)0;,"^(x,r')Vr(^,^') 



(2.6) 



where ai and a2 are replica indices, and — s- is the 
number of replicas. Of course, all other terms in the 
action also are replicated TV times. 



B. Composite variables 

We now proceed by rewriting our model in terms of 
variables that are more suitable for our purposes than the 
basic fermionic field. First of all, we decouple the spin- 
triplet interaction by means of a Hubbard-Stratonovich 
transformation. All other terms^we rewrite in terms 
of bosonic matrix fields Q and A. The latter proce- 
dure exactly follows Ref. 13, and we refer the reader to 
that reference for details. Here we just mention that A 
serves as a Lagrange multiplier whose physical interpre- 
tation is a self energy, while Q is isomorphic to bilinear 
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products of fermion fields. We perform a Fourier trans- 
form from imaginary time r to Matsubara frequencies 
cj„ = 27rr(n+ 1/2), 

JO 

Vin,a(x) = VT/ dr e-'"""Via(x,r) . (2.7a) 
Jo 

For later reference, we also define a spatial Fourier trans- 
form 



V'n,Q(k) = j dx 



^n,a(k) 



dji e'^'"" IpnA^) : 



(2.7b) 



and analogously for other position depeiidend quantities. 
The isomorphism then takes the form,E3 



Qi 



/ -i>iiip2^ -V'iTV;2i -V'iTV'2i V'iT'02T \ 

-"01X^^21 -t^ltn -V'U'/'2i '0Ul/'2T 
V^UV:2T '0U'02i lhl1p2l -tplltp2-\ 

\ -'0itV'2T -V'iTV'2i -V'iTV'2i V'1T'02T / 

(2 



Here all fields are understood to be taken at position 
X, and 1 = (ni,ai), etc., comprises both frequency and 
replica labels. It is convenient to expand the 4x4 matrix 
in Eq. ( ^.§| ) in a spin-quaternion basis. 



Qi2(x) = ^ (t^ (g) s,) ;Qi2(x) 



(2.9) 



and analogously for A. Here tq — sq — I2 is the 2x2 unit 
matrix, and Tj = —Sj = —i<7j, {j = 1,2,3), with cri,2.3 
the Pauli matrices. In this basis, i — and i = 1,2,3 de- 
scribe the spin-singlet and the spin-triplet, respectively. 
An explicit calculation reveals that r = 0, 3 corresponds 
to the particle- hole channel (i.e., products ip'4')j while 
r = 1^2_describes the particle-particle channel (i.e., prod- 
ucts tpip or ipip). For our purposes the latter will not be 
of importance, and we therefore drop r = 1,2 from the 
spin-quaternion basis. 

Denoting the Hubbard- Stratonovich field by M, we 
now can exactly rewrite the partition function as an in- 
tegral over the three fields Q, A, and M, 



Z= / D[Q,A,M] 



^A[Q,A,M] 



with an action 



(2.10a) 



-iA 



A[Q, A, M] = ^dis[Q] + A^nt [Q] + \ Tr In (g^^ 

+Tr (aq) -|dx^^^^M:(x)^M^„(x) 



Q n i—\ 



r^0,3 



^tr [(r.^s.) gjJ.VnW] • (2-lOb) 



Here Tr denotes a trace over all degrees of freedom, 
including the continuous real space position, while tr 
is a trace over all discrete degrees of freedom that are 
not summed over explicitly. The first two terms in Eq. 
( p.lObj ) read explicitly. 



-4dis[(9] = 



1 



27riVFTo: 



dx tr (Q(x))' 



(2.10c) 



A 



2np(s) 



/dx^(-i)^ E E 



ni,n2,m a 



X [tr ((r.^So)Q"r,n,+™W)] 

X [tr {{rr^So)QZ+m,nA^))] 



(2.10d) 



They are simply the last two terms in Eq. (2.5a), rewrit- 
ten in terms of the Q. Finally, 



Go" 



-dr + V^/2m + n 



(2.10e) 



is the inverse Green operator. For later reference we also 
note that if A in the Tr In term in Eq. ( ^.101 ) is re- 
placed by its saddle-point value, one obtains the inverse 
saddle-poipt Green function Ggp as the argument of the 
logarithm.ll^ For our purposes it suffices to treat the dis- 
order contribution to the self energy in Born approxima- 
tion, and to neglect the (Hartree-Fock) interaction con- 



tribution. Gsp then reads 



Gsp(k, uj„ 



2m 



2t, 



■ sgn UJ„ 



(2.10f) 



Notice that Eqs. ( |2.10 ) can also be obtained from Eqs. 
(2.22) - (2.25) in Ref. |l| by subjecting the spin-triplet 
interaction term to a Gaussian transformation, and by 
using the symmetry properties of the Q-matrices. How- 
ever, the above derivation makes it clear that the order 
parameter field M is introduced in a standard way, and 
the only difference to a standard treatment is that the 
fermionic parts of the action have been rewritten in terms 
of Q and A. 



C. Separation of soft and massive modes 

The reason for our rewriting of the fermionic part of 
the action in terms of bosonic matrix fields in the previ- 
ous subsection was that this formulation is particularly 
well suited for a separation of soft and massive modes. 



a n i—1 
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For the purpose of analyzin g a Fermi liquid, this separa- 
tion was carried out in Ref. 13, and we briefly recall the 
most important points of that prodecureJla One first uses 
group theory arguments to show that the most general 
Q can be written in the form 



Q = SPS- 



(2.11) 



Here P is block-diagonal in Matsubara frequency space, 
P={ \. 1 , (2.12) 



with 



P> and P< matrices with elements P„„i where 
n, TO > and n, to < 0, respectively. For a sys- 
tem with N replicas and n Matsubara frequencies, 
the matrices S are elements of the homogeneous space 
USp(8iVn,C)/USp(47Vn,C) x USp(4iVn, C).N As such 
they can be expressed in terms of matrices q whose el- 
ements qnm are restricted to frequency labels n > 0, 
TO < 0, 



S = 



Vl - bbt b 

-6t Vl - b'fb 



where 



with 



b{q,q^) 



1/2 



(2.13a) 



(2.13b) 



(2.13c) 



Ward identities ensure that the P are massive while the 
q are massless. The latter are the diffusive modes or 'dif- 
fusons' that we mentioned in the Introduction. Similarly, 
if A is transformed according to 



A(x) =5-i(x)A(x)5(x) 



(2.14) 



A can be shown to be massive. This representation thus 
achieves the desired separation of modes. In Ref. |l^ the 
q were the only soft modes. In the present case, the order 
parameter field M is massive in the paramagnetic phase, 
but becomes soft at criticality, so it must be handled 
together with the diffusons. 

The next step is to expand the massive modes about 
their expectation values. 



P = {P) + AP , A = (A) + AA 



(2.15) 



As explained in Ref. |3|, the expectation values (P) and 
(A) can be replaced by simple saddle-point approxima- 
tions. The saddle point used in this reference is also a 
saddle point of our current action, Eq. ( ^.10b| ), if it is 
supplemented by the saddle-point value of M, (M) — 0. 
If P and A are integrated out in saddle-point approxima- 
tion, i.e. if one neglects all fluctuations of these fields, 
then one obtains for the fermionic degrees of freedom the 



nonlinear siftma model (NLcrM) that was first proposed 
by WegneiEj as an effective field theory for the disor- 
dered electron pffjohlem, and later studied extensively by 
him and others.BO Following the same procedure here, 
we get the NLcrM, but without the spin-triplet interac- 
tion term, from the first four terms on the r.h.s. of Eq. 
( ^.101: ) , and a coupling between the order parameter field 
and both the soft and the massive fermionic modes from 
the last one. We will also need the corrections to the 
NLctM, and thus rewrite the action in the form 

A[q, A/, AP, AA] = ^NL.M [q] + 5A[AP, AA, q] 
- I dx^^^M<(x)'M!„(x) 

a n i—1 

f dx^^^^M:(x) iv^r 

J an i=i 1 

X^tr {Tr®S,) [QZ, 



(x) 



^ (5AP5-1)"" ^ (x) 



(2.16) 



Here Kt = 'nN'^Tt/2. .Anlo-m is the action of the nonlin- 
ear sigma model, 

-^NLcrM 

+2H I dx tr (nQ{x 



^ /dxtr (vQ(x)^' 



-Ai:l[7rNT.Q/2] 



(2.17a) 



with A[^^ from Eq. (|2.10d|), and 



qq' 



q^q 



(2.17b) 



in terms of q, and Q a frequency matrix with elements 
^^12 = (to «> So) Si2 w„i , (2.17c) 

The coupling constants G and H are proportional to the 
inverse conductivity, G oc 1/ct, and the specifig.hpat co- 
efficient, H (X "f = lim^^o Cv/T, respectively.B'EJtJ 

SA contains the corrections to the nonlinear sigma 
model that were given in Ref. |l^. We list explicitly the 
terms that are bilinear in the massive fluctuations AP 
and AA, but do not contain couplings between the mas- 
sive modes and q, 

6A^^^ = ^dis[AP] + y dx tr (AA(x) AP(x)) 

-f i I dxdy tr (G(x - y) AA(y) G(y - x) AA(x)) . 

(2.18) 
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D. Effective field tlieory for tlie soft modes 



x(k,r]„) = -7VF + o(k2,r!„) 



(2.20a) 



So far we have exactly rewritten the microscopic action 
in a form that separates the soft modes from the massive 
ones. In this subsection we approximately integrate out 
the massive modes to arrive at an effective action that 
is capable of describing the critical behavior at the fer- 
romagnetic transition. We will also add some terms to 
the bare action that we will find later to be generated 
by the renormalization group. 



Ill 



we will derive 



In Sec. 

these terms, and will also justify our approximations and 
show that they do not influence the asymptotic critical 
behavior. 



1. Integrating out the massive modes 

We now need to dispose of the massive modes. Clearly, 
we cannot just ig nore the mas sive fluctuations. It is ob- 
vious from Eqs. ( ^.16 ), ( 2.17 ) that they are needed to 
bring the purely M-dependent part of the action in a 
standard Landau-Ginzburg- Wilson (LGW) form. Let us 
first integrate out AP and AA in a Gaussian approxima- 
tion, while neglecting the coupling between q and these 
massive fluctuations. We will consider the effects of this 

and the last term in Eq. 



coupl ing later. From Eq. ( 2.1^ 
( ^.16D with S = 1 we obtain a contribution to the action 
that is quadratic in the order parameter field M . Com- 
bining it with the AI^ term in Eq. (2.16) yields a term 
proportional to 



EEEE^<W [i + 2ra(k,i7„)]^M!„(-k) , 

(2.19a) 



k Q n i—1 



where 

^ni— n2,n 2?nin2 (k) , 

"1 ,«2 

(2.19b) 

is a restricted spin susceptibility that is given in terms of 

1 



27r7VFrei 



(k) 



I?„m(k) = (/7„m(k) 1 

with 

</'nm(k) ^ -^^G'sp(p,W„)Gsp(p-|-k, 



(2.19c) 



(2.19d) 



Here Gsp is the saddle-point Green functio n from Eq. 
( ^.10^ . The Theta-function in Eq. ( |2.19bD , which re- 
stricts the frequency sum to frequencies that both have 
the same sign, results in x being the non-hydrodynamic 
part of the spin susceptibility of noninteracting disor- 
dered electrons. For small frequencies and wavenumbers, 
it reads 



For later reference, we also note that for nm < 0, I?„m 
is the basic diffusion propagator. In the limit of small 
frequencies and wavenumbers one finds 



2^riT?j (k) — 



27riVF 



i?k2 + |17„_,„| 

EE 2TTNFGHVn-m{k) , 



(nm < 0) 



with 



■Dnik) - 



1 



k2 + GHlnr 



(2.20b) 



(2.20c) 



Here D = 1/GH is a bare diffusion coefficient, and fin = 
2-KTn is a bosonic Matsubara frequency. 

Within the approximation, Eq. (2.20a), the effective 
action has the form 



Aff = -5IEEE'^'^»W [t + ak2 + 6|r!„|] 

k Q n i—1 

x''Af_„(-k)+^NL^M[(9] 



xtr 



a n i—1 
{Tr (8l Si) Qm"„i+„(x) 



r=0,3 



(2.21) 



Here t — 1 — 2iVFrt, a and b are constants, and in the 
first term we have truncated the gradient expansion after 
the leading wavenumber and frequency dependent terms. 
Notice that the first term taken in isolation describes a 
ferromagnetic transition at t = 0, which represents the 
Stoner criterion. 

The following structure emerges. Our effective action 
takes the form of a LGW action for the order parameter 
field M , a nonlinear sigma model for the soft fermionic 
modes, and a coupling between the order parameter and 
the composite fermion field Q. In our current approxima- 
tion, only the Gaussian part of the LGW action appears. 
However, keeping terms of higher order in AP would 
clearly produce terms of higher order in M, with coeffi- 
cients that are numbers in the limit of zero wavenumbers 
and frequencies. For our purposes it will suffice to keep 
the term of order M"*. We thus should add to the action 
a term 

"^LGW =U4 fd^Tj2 E (M"(x,c.„J • M"(x, ..„,)) 

ni, 712,^3 a 

X (M"(x, w„3) • M"(x, -UJn, - CJ„, - CJ„3)) , 

(2.22) 

with U4 a number, and M the 3- vector whose components 
are the *M. A detailed derivation shows that there is 
also a term of third order in M, with a coefficient that 
is proportional to either a gradient or a time derivative. 
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This term is irrelevant for the critical behavior at the 
ferromagnetic transition,u and we neglect it. 

Both from simple physical considerations, and from 
Ref. ^ it is obvious that another term of 0{M'^) must 
exist. This is the 'random mass' term that reflects spatial 
fluctuations in the location of the critical point, and it 
has the structure 



4gw -v^Jd^Y.!! |M"(x,c.„J| = 



ri\,n2 a, 13 

X \M^{x,LUn,)\ 



(2.23) 




with Vi SL number. To see how this term arise s in our 
present formulation, wc go back to Eq. ( 2.16| ). If we 
expand the iS in the last rtcpp in powers of q, the two 
lowest order contributionsEj'Efl have the structures 

Vf J dx Af (x) AP(x) , (2.24a) 

Vt J dx M{x) AP(x) q{x) q{x) . (2.24b) 

Here we have dropped frequency labels and all other de- 
grees of freedom that are irrelevant for power counting 
purposes. Let us depict the fields M, AP, and q by 
dashed, wavy, an d solid lines , respectively. Then the two 
vertices in Eqs. ( 2.24a , 2!^24tj ) have the structures shown 
in Fig. 0. Contracting the wavy lines yields an effective 



cr\j\r 



FIG. 1. M-AP, and M-AP-q^ vertices. Broken lines de- 
note M, solid lines q, and wiggly lines AP. 

vertex of the form 

T J dx M'^{x)q^{x) , (2.24c) 

which is shown in Fig. 0. This vertex can be used to con- 



FIG. 2. Effective M^-q^ vertex obtained by contracting the 
AP fluctuations of the two vertices shown in Fig. ^. 

struct the one-loop contribution to the term of O(AP') 
that is sho wn in Fig. |^. Calculating the diagram yields 
Eq. (2.23). We will further discus s this t erm from a 
power-counting point of view in Sec. |IIIC2| below. 





FIG. 3. One-loop renormalization of the vertex by 
means of the vertex in Fig. H. 



2. Effective action, and leading corrections 

We now can assemble our effective action by collect- 
ing the various terms derived in the previous subsection. 
Before doing so, however, it is illustrative and convenient 
to add a term to the LGW part_of the action that will 

As we will see, for 



be generated by the RG in Sec. Ill 



dimensions 2 < c? < 4 the gradient squared term in the 
flrst contribution to A is not the leading wavenumber 
dependence. Rather, as one would expect from Ref. 
there is a term proportional to [kj''"^, which first ap- 
pears at one-loop order. We therefore add this term right 
away. Notice that the resulting action is still local in the 
sense of Ref. ^. It will also turn out that the term in 
the LGW action that is linear in frequency is not the 
leading frequency dependence. Rather, the coupling be- 
tween M and Q effectively produces a term proportional 
to |17„|/k^E Physically, this term reflects the fact that 
spin is a conserved variable, i.e. the characteristic fre- 
quency scale vanishes in the long-wavelength limit even 
away from the critical point. We can therefore drop the 



frequency dependence from the first term in Eq. (2.21 



Taking all of these points into account, we obtain the 
following effective action, 

AeS = AlGw[M] + ^NLaMM + A[M, q] ■ (2.25a) 
Here .4lgw is the modified LGW part of the action, 

3 



A 



LGW 



[M] 



k a n i—1 

+a2k2] ^A^J-k) 
+^[gwM+4g'^M (2.25b) 



with Qrf_2 and 2 con stants, and -^l'q^ and -Ayt^^s, from 



^LGW 



Eqs. ( 2.22| ) and ( ^.23| ), respectively. The nonlinear sigma 
inodcl part of the action, AisihaM , has been given in Eq. 
(^.17a ), and represents the coupling between M and 

A,[M,q] = /dx^^^^Af:(x) ^ {V^y 

an 1=1 r=0,3 

X ^ tr [(r, ^ Si) Q^%+„(x)] . (2.25c) 



It is useful to deflne a field 
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6i2(x) = ^(Tr «) Si) ^^laW ' 



(2.26a) 



with components 
;6i2(k) -<5„,„,(-)'-/2^5„.„,_„, fM:^(k) 

n 

+ {^Y+''MZi^)] ■ (2.26b) 

In terms of b, the couphng part of the action can be 
written 



the theory. It is obvious, and easily confirmed by keep- 
ing a source term for the electron spin-density, that the 
expectation value of the order parameter field M deter- 
mines the magnetization m. Specifically, 



m = /iBv/2TVf;C<=o(x)) 



(2.29a) 



The two-point M-vertex is therefore proportional to the 
magnetic susceptibility. 

Other relevant observables are the specific heat coeffi- 
cient 



47r 

7 = yH , 



Ac[b,q] = ^/TrTKt j dx tr (^fe(x)(3(x)^ . (2.27) and the electrical conductivity 



We stress that we have added the |k|'' ^ term in Eq. 
( p5b| ) for convenience only. One could equally well work 
with the theory given by our action with ad-2 — 0, but 
some effects that are included in the bare ^eff would then 
appear only at loop level. We also stress that the coeffi- 
cients of all non-Gaussian terms in ^off are finite in the 
limit of zero wavenumbers and frequencies, so the theory 
is local in the sense of Ref. ^, in contrast to the situa- 
tion in Ref. ^ More importantly, the physically relevant 
point is that the current formulation, as opposed to the 
one in Ref. ^, makes obvious the existence of two time 
scales, as we will see in Sec. Ill below. 



In order to justify the omission of the terms left out of 
our effective action, we also need to know the structure 
of the omitted terms. We therefore list these here, using 
notation that leaves out anything that is not needed for 
power counting purposes. 

First of all we have the corrections to the nonlinear 
sigma model. In addition to Eq. ( 2. IS), th ey consist of 



those terms in the Tr In term in Eq. ( 2.10b| ) that contain 



cubic or higher orders of AA, the part of A\^^ that con- 
tains the massive fluctuation AP, and corre ction s to the 
saddle-point approximation for (P) in Eq. (2.15). These 
have all been discussed in Ref. and the same discus- 
sion applies here. 

Second, there are the terms that couple M, AP, and 
q, se e the second term in the bracket in the last line of 
Eq. (2.16). In general, they have the structure 



dnVT j dx M(x) AP(x) g" (x) , 



(2.28) 



with coupling constants dn and n = 0, 2, 3, . . .S The first 
two t erms in this expansion have already been given in 
Eqs. d^H). 



3. Observables 

For a physical interpretation of any results obtained 
from our effective action we need to identify the appro- 
priate observables in terms of the coupling constants of 



(2.29b) 



(2.29c) 



(see Sec. |IC|). Also of interest is the electronic density 
of states per spinll^ 



.0 , (2.29d) 

where the energy or frequency lu is measured from the 
Fermi energy ep- 



III. RENORMALIZATION GROUP ANALYSIS 

In this section, we first consider low-order perturba- 
tion theory to see how those terms in Eqs. ( 2.25 ) that 
were not in the bare action are generated. We then do a 
power-counting analysis to determine the minimal effec- 
tive action that needs to be analyzed in order to find the 
critical behavior at the ferromagnetic transition. Finally, 
we show that the terms that were omitted from the effec- 
tive action are irrelevant, in the RG sense, for the critical 
behavior. 



A. Perturbation theory 

We first set up a standard perturbative expansion for 
our effective action, starting with the Gaussian theory. 



1. Gaussian propagators 



We start by expanding the effective action, Eqs. ( 2.25 ), 
to bilinear order in M (or b) and q. We obtain for the 
Gaussian action 

Ao[M, 9] = - E E E E "2(k) ^M!„(-k) 



k n a i—1 



-?^E E E''?12(k)^rg34(k):.<Z34("k) 



G 



k 1,2,3,4 i,r 



-4y^EEE^'?i2W>i2(-k) , 

(3.1a) 



k 12 i,r 
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where 

U2(k) = t + ad-2m'^~^ + a2k'^ 
The bare two-point g-vertex reads 

^ ■ - ■ ■ ■ ■ -n2 j 
-<^l-2,3-4<^QiQ2'^aiQ3 2nTGKs 



0p(2) 



12,34(k) = '^13'524 (k^ + GHfln,-r 



l,2,3p(2) 



12,34 (k) — '^13<524 (k + GH^m —712 ] 



(3.1b) 



(3.1c) 
(3.1d) 



with Ks = -7riV|r,/2. 

The quadratic form defined by this Gaussian action is 
easily inverted. For the order parameter correlations we 
find 

(^M"(k) Wl(p)) = ^k,-p(5„._™5,, S^p^Mnik) , 

(3.2a) 

(r^l2(^) i^34(P)) ^ p ^rs ^ij ^aia2 ^cxicts 

xMni-n2ik) [(^1-2,3-4 - (-)''<5l-2.4-3] , 

(3.2b) 

in terms of the paramagnon propagator 

1 



t + a,_2|k|''-2 + a2k2 + ^.';^-g^l I 



(3.2c) 



Notice that the coupling between the order parameter 
field and the fermionic degrees of freedom has produced 
the dynamical piece of M that is characteristic of disor- 
dered itinerant ferromagncts.0 

For the fermionic propagators we find 

(r9l2(k) ig34(P)> ^ K-P ^rs 5i] ^ *rf2^_3/(k) , (3.3a) 



in terms of the inverse of Eq. ( 3.1c ) 



^ 12,34 



x27rTGi^,I?„,_„,(k)I?i^^L„^(k) , (3.3b) 

and 

' ' ^12,34 (k) — di3 524 T^ni —7X2 (k)-5i-2,3-4 5 0.10.2 *^Ql Q3 

x2T:TGKt (!?„,_„, (k))' X„,_„,(k) . 

(3.3c) 

Here P^^^ is the spin-singlet propagator, which in |the 
limit of long wavelengths and small frequencies reada3 



P(f)(k) = 



k2 + G{H + Ks)n.r, 



(3.3d) 



Finally, due to the coupling between M and q we have 
a mixed propagator 



f ^rs ^ij ^aia2 ^aic 



G 



(k)X„,_„,(k) [(5l-2,3-4 + (-)''+'<5l-2,4-3] ■ 

(3.4) 




FIG. 4. One-loop renormalization of the vertex U2. 



2. One-loop order 



Let us now consider the renormalization of the Af^- 



vertex U2, Eq. ( 3.1b ). At one- loop order, the relevant 
diagram is the one shown in Fig. |^. While the complete 
result is rather involved, it can be simplified by means of 
the following observation. The structure of the diagram 
leads to a frequency-momentum integral over diffusion 
poles multiplied by one or more paramagnon propaga- 
tors. Inspection shows that the frequency in this integral 
scales like a wavenumber squared. To leading order in 
the distance from the critical point, and for the purpose 
of obtaining leading infrared singularities, we therefore 
can use the following approximation in the integrand, 



GKt^i P,(p)A4,(p) = l + 0(t,|p|'') 



(3.5a) 



At zero external frequency the diagram then yields 



5u2(k) 



2 V 



^27rT^I?;(p)I?;(p + k) . (3.5b) 



1=0 



Performing the integrals shows that this diagram pro- 
vides a finite renormalization of the coefficient ad--2 of 
the |k|''~^ term. In particular, it would generate this 
term if it were not present in the bare action. This is the 



reason why we added this term by hand in Sec. II D 2 



We now consider the one-loop renormalizations of the 
other vertices in the effective action. Of particular in- 
terest are the coupling constants H and G in the two- 
point g-vertex, which determine the specific heat coeffi- 
cient, and the conductivity, respectively, see Sec. [ID 3 
above. The topological structure of the relevant diagrams 
is shown in Fig. |^. From extensive work on the nonlin- 




(a) 




FIG. 5. One-loop diagrams that renormalize H and G. 

ear sigma model it is known that diagram (b), combined 
with suitable contributions from diagram (a), is finite for 
d > 2J3 A calculation of diagrams (c) and the remaining 
parts of (a) reveals that their contribution to the one- 
loop renormalization of G consists of two pieces which, 
at zero external frequency, are given by 



(^G) 1 = ^ ^ 5: 27rT 5: P,^ (p) M (P) 

P 1=0 



(3.6a) 
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and 



^ OO 



k=0 • p 

X [Vi{p)-Vi{p + k)] . (3.6b) 

respectively. Both of these integrals are finite in d > 2. 
A simple calculation shows that the one-loop correction 
to the density of states is given by the same integral as 
((5G)i. For the iJ-renormalization, diagram (b) is again 
finite, while the other two contributions yield 

SH =lGKt:^27TTj2^T.^iip)Miip) , (3.7a) 



V 

1=0 p 



where n is the external frequency label. For later refer- 
ence we note that diagrams (a) and (c) in Fig. H each 
contribute one half of this result. Individually, each of 
these diagrams also contributes pieces that diverge like 
l/fi„; these contributions cancel between the two dia- 
grams. (The same is true for diagram (b) and the cor- 
responding contributions from diagram (a).) Notice that 
the frequency structure is slightly different than in the 
case of the G-rcnorm aliza tion. This leads to a finite fre- 
quency sum in Eq. ( ^.7a ), and as a result the integral 
is logarithmically infrared divergent for all dimensions 
2 < d < 4, 



SH=-GKt ln(l/r!„) 



(3.7b) 



Here G — G Sd/i^ir)'^, with Sd the surface area of the 
unit {d — l)-sphere. While this is consistent with the 
result of Ref. | that the specific heat coefficient at the 
quantum ferromagnetic transition is logarithmically di- 
vergent for all 2 < d < 4, it is obvious from the present 
formulation of the theory that this is very unlikely to 
be the exact critical behavior as claimed in that refer- 
ence. The reason is that an insertion of the one-loop H 
into the one-loop diagrams yields a In'^ Q singularity, and 
so on. Similarly, inserting the one-loop result for the q- 
propagator in the diagram shown in Fig. ^ generates a 
logarithmic correction to the vertex U2(k). Unless these 
insertions are exactly cancelled by skeleton diagram con- 
tributions, the exact critical behavior must therefore in- 
volve a more complicated function of lnf2. We will come 
back to this point in the next subsection, and in II. 

Finally, the diagrams shown in Fig. |^ also determine 
the renormalization of Ks- Again, diagram (b) yields 
a finite result, while from the other two diagrams one 
expects 



SK, = -SH 



(3.7c) 



with 6H the divergent part of the i7-renormalization, Eq. 
( |3.7b| ). This result follows from what is known about the 



nonlinear sigma model as we will discuss in Sec. IV A. 

Notice that the effects we have discussed above at one- 
loop level would of course also occur if we had worked 



with the theory that one obtains if one puts the bare 
O'd-2 equal to zero. However, in that theory their deriva- 
tion would have required a renormalization of the para- 
magnon propagator. In a diagrammatic language, by 
adding the Ikl**^^ term to our bare action, we have made 
use of skeleton diagrams that contain certain infinite re- 
summations. We note in passing that one might worry 
about higher orders in the loop expansion producing even 
stronger nonanalyticities than jk]'^^^. We will show in II 
that this is not the case. 



B. Naive fixed points and their instability 

We now proceed to perform a power counting analysis 
of our effective action, Eqs. (2.25). Our goal is to under- 
stand the perturbative results of the preceding subsection 
from a more general point of view, and to determine the 
minimal effective action which, when solved, will yield 
the exact critical behavior. For this purpose it is conve- 
nient to rewrite the action in a schematic form which sup- 
presses everything that is not necessary for power count- 
ing, 



+0[diM^,M^) 



G 



dx{d^qy + H J dxVLq^+KsT J dxq^ 

j dxdU^ + Hi J dxilq* 
+0{Tq\dlq^^lq'') 
+VTciJ dxMq + VTc2 J dxMq^ 

+0{VfMq^) . (3.8) 



Here the fields are understood to be functions of position 
and frequency, and only quantities that carry a scale di- 
mension are shown. The bare values of G4 and Hi are 
proportional to those of G and iJ.E3 The term of order 
Tq^, which arises from the interacting part of -Anlo-m, 
will not be of importance for our purposes although its 
coupling constant squared has the same scale dimension 
as I/G4 and H4. It is therefore not shown explicitly. 
The ci, C2, etc. are t he cou pling constants of the terms 
contained in Ac, Eq. ( ^.25c ). 

We now assign a length L a scale dimension [L] = 
— 1. Under a renormalization group transformation that 
involves a length rescaling by a factor 6, all quantities 
will then change according to A ^ b^^^ A, with [A] the 
scale dimension of A. In particular, imaginary time t 
and temperature T or frequency Q have scale dimensions 
[t] = — z, and [T] = [Q] = z, respectively. 



10 



1. Hertz's fixed point 

To illustrate an important point, let us first show how 
one recovers Hertz's mean-field fixed pointla within the 
present formalism. Let us look for a fixed point where 
the coefficieiits 02 and Ci are marginal, [02] = [ci] ~ 0. 
This choicec3 is motivated by the desire to find a fixed 
point with mean field-like static critical behavior, and 
with dynamics given by the frequency dependence of the 
standard paramagnon propagator , Eq. (3.2c). (Recall 
that the frequency dependence of M. was produced by the 
vertex with coupling constant ci.) From the condition 
that the action be dimensionless we then obtain the scale 
dimensions of the order parameter field, 

[M„(x)] - (d-2)/2 , (3.9a) 

and wc find t to be relevant with scale dimension 



M=2 



(3.9b) 



The correlations of the q-field we expect to descpibe the 
diffusive dynamics of the fermions, so we choosdlj 



[g„™(x)] = (d-2)/2 



(3.9c) 



and G, i?, and are all dimensionless. The marginality 
of the coupling constant ci then implies 



z - [T] = 4 



(3.9d) 



This is the fixed point proposed by Hertz,! which leads 
to mean-field critical behavior. It is unstable because the 
coupling ad-2 is relevant with respect to this fixed point, 
as has been pointed out in Ref. ^. While this is obvious 
from the action as formulated here, the following inter- 
esting question arises. Suppose we had not added the 
term with coupling constant ad-2 to our bare action. 
Since this term was generated by means of the Mq^- 
vertex with coupling constant C2, see Fig. ^, C2 should 
be relevant with respect to Hertz's fixed point. However, 
power counting with the above scale dimensions yields 



[C2] = ^ (d + z - 6) 



(3.9e) 



so with the above value z = 4, C2 seems to be irrelevant. 
The resolution of this paradox lies in the fact that there 
is more than one time scale in the problem, and hence 
all factors of T do not carry the same scale dimension 
z. This is obvious if we consider the fermionic sector 
of our action; the factors of T in the q^-vertex carry a 
scale dimension [T] — Zdig = 2, which corresponds to 
the diffusive time scale that describes the dynamics of 
the electronic soft modes and which is distinct from the 
critical time scale that corresponds to \T] = Zc = 4. The 
scale dimensions of the factors of VT in the coupling 
part of the action are therefore not a priori clear, and 
they may depend on the diagrammatic context a vertex 
is used in. Consider the diagram in Fig. || again. In 



this context, the two factors of -s/T contribute to the 
frequency measure of a fermionic loop, and hence they 
carry the diffusive time scale. Indeed, with z = 2 Eq. 
( |3.9(j ) shows that C2 is relevant for 2 < < 4, and its 
scale dimension is consistent with that of 

It is worthwhile to mention that o ne ca n tell a priori 
that the scale dimension for i, Eq. ( 3.9b| ), cannot cor- 
respond to a stable fixed point since it corresponds to 
a correlation length exponent v — l/[t\ — l/2,pjjvhich 
violates the Harris criterion inequality v > 2/dLJ The 
relevance of C2 provides an explicit mechanism for the 
instability. 

We also note that the above discussion is oversimplified 
in that it pretends that M has always the same scale di- 
mension, independent of the context the order parameter 
fiuctuations appear in. As we will see in the next sub- 
section, this is not quite true. However, since this point 
is not crucial for the instability of Hertz's fixed point we 
have suppressed it. 



2. A marginally unstable fixed point 

Given the presence of the term with coupling constant 
ad-2, an obvious attempt to find a stable fixed point is 
to jdjioose ad-2 and ci to be marginal instead of 02 and 
ci.Ea A slight complication, however, lies in the fact that 
due to the existence of two time scales, ad-2 will not nec- 
essarily be marginal under all circumstances. Na mely, if 
the frequency in the paramagnon propagator, Eq. (3.2c), 
is diffusive, i.e. if it scales like k^, then ad-2 will be irrel- 
evant. As we will see below, this can happen if the para- 
magnon propagator appears as an internal propagator in 
perturbation theory, although in the critical paramagnon 
ad-2 is marginal. In general, we therefore demand only 
that Ci be marginal, that the scale dimension of the q- 
field be consistent with a diffusive (qq) propagator. 



['?r™(x)] = i(d-2) , 



(3.10a) 



and that the diffusive time scale be represented by a dy- 
namical critical exponent 



Equation (3 



^^diff — 2 

then implies 

G] = [H] = [Ks 







(3.10b) 



(3.10c) 



The marginality of ci implies for the scale dimension 
of the order parameter field 



[M„(x)] = l + (d-z)/2 



(3.11a) 



with z the dynamical exponent associated with the %/T 
prefactor in the ci-vertex. In the critical paramagnon 
propagator we expect ad-2 to be marginal, which implies 
[Af„(x)] = 1, and hence a critical time scale characterized 
by 
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Zr = d 



(3.11b) 



and a critical exponent rj, defined by [M,i(x)] = (d — 2 



T] = A - d 



(3.11c) 



This makes 02 irrelevant, while t is relevant with scale 
dimension 



[t]=d-2 , 



(3.12a) 



which leads to a correlation length exponent 

V =l/[t] = l/{d~2) . (3.12b) 

Notice that in contrast to the situation at Hertz's fixed 
point, this result respects the Harris criterion.EZlHere and 
in the remainder of this paper we restrict ourselves to 
the range of dimensions 2 < d < 4, which includes the 
physically interesting case d = 3. For the behavior in 
higher dimensions, see Ref. |^. 

The preceding results characterize the Gaussian fixed 
point that was discussed in Ref. ^ If all other terms in 
the action were irrelevant, or marginal leading to finite 
renormalizations only, then this fixed point would be sta- 
ble. To check this, we need to consider the corrections 
to the Gaussian action. We start with C2, whose scale 
dimension is 



[C2] = 1 - z/2 



(3.13) 



with z the scale dimension of the factor of T in that 
vertex. If this temperature represents the critical time 
scale, then C2 is irrelevant. However, if it represents the 
diffusive times scale, then it is marginal. This can indeed 
happen, as we have discussed in Sec. IIIB 1 above. The 
example we used, viz. the diagram in Fig. |, just leads 
to a finite renormalization of the coefficient ad-2i which 
is part of our effective action anyway. If this were the 
only effect of C2, then we could neglect it. However, this 
is not the case. The o ne-loop renormalization of H that 
was discussed in Sec. [II A 2 provides an example how 
operators that appear irrelevant by naive power counting 
can be effectively marginal due to the existence of two 
time scales, lead to logarithms, and therefore need to be 
kept. Consequently, C2 is not necessarily harmless even 
if z = d in Eq. (|3T|). This is an important point which 
we now discuss in detail. 

Consider diagrams (a) and (c) in Fig. ||. They both 
lead to a correction to the two-point g-vertex that is of 
the form, at zero external wavenumber and frequency. 



5T 



(2) 



^-^E ^E^'(P)-^'(P) ■ (3.14a) 



1=1 



For scaling purposes, let us cut off the momentum in- 
tegral in the infrared by 1/6, with h a RG length scale 
factor. Doing the integral then shows that it is given 



by a constant plus a term proportional to 6 In b. The 
constants cancel between the two diagrams, and we have 



In 5 



(3.14b) 



Notice that the frequency in the above integral scales like 
a wavenumber to the power d, so (02)^ in diagram (c) 
has a negative scale dimension — (d — 2), and so does the 
quartic vertex in diagram (a). The salient point is now 
as follows. For the purpose of the renormalization of G, 
i.e. the wavenumber dependent part of T^^\ we need to 
replace 1 /b by k. We then obtain the gradient squared of 
the bare vertex times a factor |k|''^^. The contribution is 
therefore irrelevant, in agreement with the negative scale 
dimensions of the vertices and the result of the explicit 
perturbative calculation. However, for the purpose of 
the renormalization of H we need to replace 1/5 by an 
appropriate power of the frequency. This can be Sl^/^, 
if represents the diffusive frequency scale, or fl^^'^, if 
it represents the critical one. Since the frequency in the 
integral scales like |p|'', the latter applies and we have 



(3.14c) 



in agreement with Eq. ( 3.7b ). 

The point illustrated above is as follows. Due to the 
existence of two different time scales, the fact that an 
operator has a negative scale dimension by naive power 
counting, which is based on the consideration of length 
scales, does not necessarily imply that it will be irrele- 
vant. Rather, operators with scale dimensions between 
zero and — (d — 2) may act as marginal operators with 
respect to frequency scaling. Notice, however, that for 
this mechanism to be operative it is crucial that the 
vertex being renormalized is proportional to frequency. 
Therefore, the seemingly irrelevant operators become ef- 
fectively marginal with respect to H , but not with respect 
to G or any other coupling constant. In the appendix we 
discuss another aspect of this phenomenon. 

We conclude that the Gaussian fixed point of Ref. ^ is 
not stable since there are operators that are effectively 
marginal with respect to it. If these operators just led 
to finite renormalizations, this would still not change the 
conclusions of the earlier paper. However, as we have 
seen above, they lead to logarithmic corrections to power- 
law scaling and hence need to be kept. The problem is 
less severe than in the case of Hertz's fixed point, how- 
ever, since now there are no relevant operators. If one can 
show that all other terms are truly irrelevant, then the 
conclusion would be that to determine the exact critical 
behavior it suffices to keep the Gaussian action plus the 
M coupling and all terms up to 0{q'^). We investigate 
this hypothesis next. 



C. Effective action for the critical behavior 
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1. The effective action 

From the discussion in the preceding subsection we in- 
fer an educated guess for an effective action that contains 
only the terms needed for a description of the critical 
fixed point and the associated critical behavior. This ac- 
tion should cont ain all of the terms that are shown explic- 
itly in Eq. (3^), except that for 2 < c? < 4 one can drop 
the gradient squared term in the LGW part of the action. 
Notice that we need to keep the terms of 0{d^ q^, VI q^) 
in the expansion of the nonlinear sigma model in powers 
of q, as they give rise to diagram (a) in Fig. ||. These 
terms appear irrelevant by naive power counting, but 
contribute to the leading frequency depend ence by means 
of the mechanism discussed in Sec. [II B 2 and in the ap- 
pendix. By the same argument one should keep the terms 
of order O(g^) and 0{q'^) that arise from the spin-singlet 
interaction. However, by themselves these vertices give 
only rise to diagrams that are finite in d > 2,Q and com- 
bined with C2 or other vertices th at c ontain M they lead 
to mixed {hq) propagators, Eq. (3.4), which are less in- 
frared diver gent then the second term on the right-hand 
side of Eq. ( [3.3c| ). These terms can therefore safely be ne- 
glected. This leaves the spin-singlet interaction constant 
Ks entering the effective theory via the vertex F*^^) only. 
Since Kg ^ Q does not change the diffusive structure of 
the noninteracting q propagator, it can be dropped there 
as well. Restoring all indices, the suggested effective ac- 
tion for describing the critical fixed point readaHj 

k,n,Q: i—1 

E E'-9i2(k)^rg34(k);<z34(-k) 

k 1,2,3,4 rA 

~4G EE E y E 

1,2,3,4 r,S,t,U 21,22,23,24 ki,k2,k3,k4 

X^^^=.:f,!r;y34(ki, k2, kg, k4) ;^<Zi2(ki) :^(732(k2) 



The bare values of the coupling constants ci and C2 are 
related, and given by 



-ci^EE>i2(k);'?i2(-k) 



k 12 



+-^^4E E EEEE^CJk) 

^ k,p ni ,"2 ,™ i',s,t 2=1 j,k a,p 

'i(p)?9^f™(-p-k)tr (t^t.t/) tr (s2Sj4) 

-ii(P) tl'^Ji'n, (-P - k) tr (r^r^Vf) tr (^s^s]sk^ , 

(3.15a) 

with F(2) from Eqs. ( |3.1cD and (|3.1(j ) with K, = 0, and 



*'T//2rr^234(ki,k2,k3,k4) - -(5ki+k2+k3+k4,0 

xtr {TrTlTtrl) tr (sj^sj^Sig^^^ (ki • kg + ki • k4 
+ki-k2+k2-k4-GiIl]„i_„J .(3.15b) 



ci = 16c2 = Ay^T: Kt 



(3.15c) 



Notice that this action is not Gaussian, and therefore the 
critical behavior is not easy to determine. 



ee will solve the effective model given by Eqs. (3.15) in 
where we will show that the exact critical behavior 
differs from the Gaussian one by logarithmic corrections 
only. In the remain der o f this paper we show that the ac- 
tion given by Eqs. (3.15) really is sufficient for describing 
the critical behavior in 2 < d < 4. 



2. Corrections to the effective action 

We now sh ow th at all terms that were neglected in 
writing Eqs. ( |3.15 ) are irrelevant by power counting, 
keeping in mind the complications due to the two time 
scales that were dicussed in Sec. [II B above. In addition 



t o the scale dimensions of M and q given in Eqs. ( 3.10a ), 
( |3.11a ), we need for this purpose the scale dimension of 
the massive fields AP and AA. The correlations of AP 
are short ranged, and of the same nature as at the Fermi 
liquid fixed point that was discussed in Ref. ^ We thus 
choose 



[AP(x)] = [AA(x)] = d/2 



(3.16) 



Power counting now proceeds as usual. The nonlinear 
sigma model action we have kept up to O(g^). Higher 
order corrections have the same scale dimensions as at 
the Fermi liquid fixed point in Ref. |l^. They thus are 
all irrelevant with scale dimensions that are smaller than 
— {d— 2) and are therefore harmless. The couplin gs be - 
tween M, q, and the massive modes given in Eq. ( 2.28 ), 
for even powers of q, have scale dimensions 



W = -^(rf-2)-|<^(d + z-2) , (3.17) 

for n > 2, and thus can safely be neglected. For odd 
powers of q, the couplings contains an effective external 
frequency, and therefore are even less irrelevant than Eq. 
( |3.17 ) suggests. Ip-particular we confirm that di, which 
we have dropped,0 has a scale dimension [di] = — 3z/2 
and is thus more irrelevant then d2- [do] = (d — 2 — z)/2, 
which becomes marginal in d = 4 if z = Zdiff = 2. How- 
ever, d = 4 is a special dimension anyway, and for d > 4 
one obtains a different fixed point since the |k|''~^ term 
in the LGW action is no longer leading. A remaining 
question is whether the formally irrelevant do can be pro- 
moted to marginal or relevant status by the same mech- 
anism that is operative for, e.g., C2- The answer is neg- 
ative, since the mechanism works only for the renormal- 
ization of H, and in order to renormalize H, dg needs 
to be combined with some d„ with n > 2. However, 
[do d2] = -z < -(d - 2) for d < 4. Therefore, aU of the 
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dn can be safely neglected. Similarly, all terms of higher 
than quadratic order in M are irrelevant. We mention, 
however, that in the ordered phase the term of 0{M'^) be- 
comes dangerously irrelevant and needs to be kept. This 
will be important in II. 



IV. DISCUSSION 



Finally, the random mass term -A^'^q'^, Eq. ( 2.23 ), de- 
serves an extra discussion. The scale dimension of the 
coupling constant in Eq. (2.23) is [U4] = d — 4, while 
the s cale dimension of the combination of dp and d2 in 
Eq. (^.28| ) which produce V4 in perturbation theory (see 
Fig. y) is [(do ^2)^] = —22: = —4. [^4] is thus much less 
irrelevant than one might expect from naive power count- 
ing. The resolution of this discrepancy is as follows. The 
?;4-vertex shown in Fig. |^ at zero wavenumber has the 
schematic structure 



dy / du;iu{q\^)q'{y)) 



(3.18) 



which has a naive scale dimension of d (with z — 2). 
However, the integral is a finite number, and so its ac- 
tual scale dimension is zero. If we consider the vertex 
function at a finite wavenumber |k| and perform a gradi- 
ent expansion, then we obtain an expansion of the form 



const. 



(3.19) 



What happens here is that power counting yields the 
scale dimension of the first nonanalytic term in the gra- 
dient expansion, but misses more dominant analytic con- 
tributions. This is of no consequence as long as the latter 
just renormalize existing terms in the action. Here, how- 
ever, they produce a new term in the action, viz. the 
random mass term, and therefore need to be taken into 
account. The difference betw een the naive scale dimen- 
sion of the integral, Eq. ( 3.18 ), viz. c?, and its actual scale 
dimension, viz. zero, is precisely the difference between 
[vi] = d - 4 and [(do ds)^] = -4. 

We finally come back t o th e simplifications inherent in 
our starting point, Eqs. (2.2), which describes the para- 
magnetic phase as a disordered electron fluid while ne- 
glecting band structure and other features of solids. The 
justification for these simplifications is as follows. The 
disordered Fermi liquid fixed point is characterized by 
relatively few parameters. E£l This is in contrast to a clean 
Fermi liquid, which requires an infinite number of Fermi 
liquid parameters, or a whole function, to completely 
characterize the fixed point .EJ The crucial physical dis- 
tinction is that for the disordered case the slowest, and 
therefore dominant, modes are diffusive and arise only 
from electron number density, spin density, and particle- 
particle density variables. In contrast, in the clean case 
there are an infinite number of soft single-particle and 
two-particle modes. This simplification for the disor- 
dered case carries over to the description of the ferro- 
magnetic quantum phase transition. 



As we have seen in Sec. Ill C , the effective action for 



the critical behavior is not Gaussian, and therefore a de- 
termination of the critical behavior is nontrivial. It turns 
out that the critical behavior in all dimensions d > 2 can 
nevertheless be determined exactly, and is given by the 
power laws found in Ref. |^ with additional logarithmic 
corrections to scaling. Thia-solution of the effective ac- 
tion will be deferred to ll£3 Here we restrict ourselves 
to a discussion of some general features of our effective 
theory, and of its relation to previous approaches to the 
problem. 



A. Relation to other approaches 

Let us briefly discuss the relation between our current 
approach and previous theori es. This is most easily done 
by starting from Eq. ( ^.10b| ). By formally integrating 
out the fermions, i.e. the fields Q and A, from this for- 
mulation of the action one obtains an LGW theory or 
action entirely in terms of the order parameter field. If 
the fermions are integratedj-put in tree approximation, 
one recovers Hertz's theory.0 If they are integrated out 
formally exactly, the vertices of the LGW functional are 
given in terms of spin-density correlation functions for 
a 'reference ensemble' or fictitious electron system that 
has no bare spin-triplet interaction. This is the theory 
that was analyzed in Ref. ^. The disadvantage of that ap- 
proach is that the reference ensemble contains soft modes, 
viz. the q, and integrating them out produces effective 
vertices in the LGW theory that diverge in the limit of 
small wavenumbers and frequencies. That is, one obtains 
a nonlocal field theory. Furthermore, Ref. ^ performed 
a power-counting analysis only, and integrating out the 
fermionic degrees of freedom obscured the subtleties that 
arise in this context due to the existence of the diffusive 
time scale in addition to the critical one. As a result, 
the power counting analysis of Ref. ^ was insensitive to 
the logarithmic corrections that we foun d by m eans of 
explicit perturbative calculations in Sec. HI A 2| and ex- 
plained in Sec. [II B 2 in terms of a more sophisticated 



scaling analysis than the pure LGW theory allowed for. 
Notice that in some other respects Ref. ^ was actually 
more sophisticated than the present theory. For instance, 
it included in the bare action effects that require a one- 
loop analysis in the present approach, e.g. the |k|'^~^ 
term in the vertex U2. However, the insensitivity to loga- 
rithmic corrections is hard to overcome within the frame- 
work of the nonlocal theory 

The relation between the present t heory and Ref. ^ is 
less obvious. To see it, consider Eqs. ( 2.25| ) and integrate 
out M . This yields a nonlinear sigma model with a triplet 
interaction amplitude that is given by the static param- 
agnon propagator. We have performed explicit calcula- 
tions within this theory, and ascertained that it yields 
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the same results as the coupled M-q theory discussed 
above, as it should. This equivalence between th e M- q 
theory and the sigma model is the basis for Eq. ( 3.7c ), 
since within the sigma model H + Kg is not singularly 
renormalized.Q The bare nonlinear sigma model in Rcf. 
^ had a point-like spin-triplet interaction amplitude, but 
under renormalization the |k|''~^ that is characteristic of 
the static paramagnon is generated. It is thus plausi- 
ble that the pure nonlinear sigma model should contain 
the critical fixed point for the ferromagnetic transition. 
However, since the order parameter has been integrated 
out, the nature of the transition is completely obscured 
within this approach, and a description of the ordered 
phase is not possible. This is the reason why Ref. ^ 
could (pJy conclude that the transition is of magnetic 
nature.E2l We will come back to the detailed connection 
between the two approaches in 11. Here we just mention 
that the present analysis positively identifies the runaway 
flow that is encountered in the nonlinear sigma model in 
the absence of any spin-flip mechanismaJ as signaling the 
ferromagnetic transition.Q 

We also mention that the flxed point identified in Sec. 



HI C above violates-some of the general scaling laws ob- 
tained by Sachdev.0 As has been discussed in Ref. ^ in 
some detail, this can be traced to the presence of dan- 
gerous irrelevant variablesp-which can always invalidate 
general scaling arguments. E3 

We conclude that all of the previous approaches to the 
problem break down at some level, and that the basic 
problem is always the same, namely a lack of explicit- 
ness. Only a local field theory that correctly identifies 
and keeps all of the soft modes allows for the explicit cal- 
culations necessary to check more general arguments that 
may break down because of the failure of hidden assump- 
tions. Interestingly, as we will show in 11, the problem 
was solved technically correctly in Ref. ^, but the miss- 
ing physical interpretation rendered this result of limited 
value at the time. 



B. Scaling issues 

Let us finally come back to the issue of the two different 
time scales, which has been crucial for a correct applica- 
tion of scal ing ideas to the problem. As we have seen in 
Sec. IIIBl, the implicit assumption of the existence of 



only one time scale, namely the critical one, can lead to 
wrong conclusions if one relies strictly on power count- 
ing arguments. Explicit loop calculations, on the other 
hand, reveal the fallaciousness of the assumption by pro- 
ducing terms in the action that are inconsistent with the 
power counting. The point is that the diffusive modes, 
whose time scale is different from the critical one, pro- 
duce long-range correlations everywhere, not just at the. 
critical point, as has been discussed in detail elsewhere.ES 
These long-range correlations are reflected, for instance, 
in the |k|''~^ term in the LGW part of the action, Eq. 



( ^.25b| ), which is responsible for the instability of Hertz's 
flxed point. 

For the instability of the Gaussian flxed point of Ref. 
^ a similar mechanism applies, although it is weaker and 
l ess o bvious. As we have seen in the context of Eqs. 
( |3.14 ), non-Gaussian terms that formally have a negative 
scale dimension can effectively become marginal with re- 
spect to frequency dependent coupling constants. This 
"counting accident" can only happen for vertices that 
vanish at zero frequency, and it has been analyzed from 
a RG point of view in Sec. HI B 2 and in the appendix. 



We will come back to these arguments in II, where we 
will provide both a resummation of perturbation theory 
to all orders and a complete scaling description of the 
exact critical behavior, including all logarithmic correc- 
tions. 
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APPENDIX A: CONSEQUENCES OF TWO 
DIFFERENT TIME SCALES 

In this appendix we discuss two additional aspects of 



the crucial point made in Sec. Ill B 2 



Let us first take a phenomcnological scaling point of 
view. The scaling equation for the two-point g-vertex 
function F^^^^ reads 

(Al) 

Here the ellipses denote the dependence of the scaling 
function 7'^' on all operators that are not shown explic- 
itly. Among these are I/G4 and H4, which play the same 
role for scaling as does. For simplicity we restrict the 
discussion to the effect of the latter. In writing Eq. (|Al| ), 
we have allowed for a dependence on both the diffusive 
and the critical frequency scale. At zero-loop order, 7^^^^ 
depends only on the former, and by putting b — 1/Vf2 
we have F(^^(k = 0, fi) cx H,. At one-loop order, it de- 
pends on the critical frequency as well, which opens the 
possibility of a stronger frequency dependence propor- 
tional to fP^'^. However, the one-loop contribution has 
the property ^^'^^ {x,y, z) = f{yz'^), which restores the 
linear frequency behavior of F*^^^. This is the same phe- 
nomenon that we have discussed within the context of ex- 



plicit perturbation theory in connection with Eqs. ( 3.14 ). 
Logarithmic corrections to scaling are neglected in this 
simple argument. 

To illustrate the same point from a RG fiow equation 
point of view, and at the same time see the origin of the 
logarithms, we absorb the freque ncy or temperature fac- 
tors multiplying H and ci in Eq. (|3.8|) into these coupling 
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constants by defining H — fffi, and Ci = ci a/T- For _ff , 
ci, and C2 we then have the flow equations 



dH 
dlnb 

dci 
din 6 

dC2 

din 6 



= 2H + const. X c^^ 

= dci/2 , 

= -(d-2)c2/2 , 



(A2a) 
(A2b) 
(A2c) 



plus higher loop orders. Again, I/G4 and play a 
role analogous to c^, and we have suppressed them for 
simplicity. The solution of this system of flow equations 



H{b) = H{b = 1) 6^ + const, x 6^ ln6 



(A3) 



In this picture, the positive scale dimensions of H and 
ci reflect the fact that frequency or temperature is a rel- 
evant variable. The critical frequency is more relevant 
than the diffusive one, but this difference is made up for 
by the fact that the critical frequency is always multi- 
plied by C2^. In this way the formally irrelevant C2 ef- 
fectively acquires a marginal status. The logarithm, at 
one-loop order, reflects a resonance between the scale di- 
mensions of ci and C2 , and represents one of the possibil- 
ities in Wegner's classification of logarithmic corrections 
to scaling|2il as was already pointed out in Ref. ^. At 
higher loop order, however, additional logarithmic terms 
appear as we will show in 11. 



^ See, e.g., S. Sachdev, Quantum Phase Transitions, (Cam- 
bridge University Press, Cambridge 1999). 

^ J. A. Hertz, Phys. Rev. B 14, 1165 (1976), and references 
therein. The finite temperature properties of this model 
were discussed by A.J. Millis, Phys. Rev. B 48, 7183 (1993). 

^ S. Sachdev, Z. Phys. B 94, 469 (1994). 

* T. R. Kirkpatriclc and D. Belitz, Phys. Rev. B 53, 14364 
(1996). 

^ T. Vojta, D. Belitz, R. Narayanan, and T.R. Kirkpatrick, 

Z. Phys. B 103, 451 (1997). 
^ For a review, see, D. Belitz and T. R. Kirkpatrick, Rev. 

Mod. Phys. 66, 261 (1994). 

In low-order perturbative RG treatments, the existence of 
this transition is signaled by runaway flow. This has led 
to widespread confusion, and even to a popular belief that 
the entire metal-insulator transition problem in the pres- 
ence of interactions is pathological. To avoid exacerbating 
this situation, we emphasize that (1) the runaway flow oc- 
curs in only one universality class, while all others show 
a perturbative fixed point that describes a metal-insulator 
transition, and (2) the technical problem disappears, and 
a critical fixed point is found, if a resummaticj«|-pf the per- 
turbation theory to all orders is performed.ElS However, 



the nature of the transition could not be identified with 
the methods employed in these papers. 
* T. R. Kirkpatrick and D. Belitz, Phys. Rev. B 45, 3187 
(1992). 

^ We call a theory "local" if the vertices in a LGW functional 
exist in the limit of zero frequencies and wavenumbers. We 
stress that this definition in general makes sense only at the 
level of a bare theory, as renormalization usually produces 
singular wavenumber dependences of sufficiently high order 
vertices. 

To avoid misunderstandings, we point out that the Gaus- 
sian theory of Refs. contains efltects whose incorporation 
within Hertz's approach requires loops. 
R. Narayanan, T. Vojta, D. Belitz, and T.R. Kirkpatrick, 
Phys. Rev. B 60, 10150 (1999). 

D. Belitz, T.R. Kir kpatrick, Maria Teresa Mercaldo, and 



Sharon L. Sessions, cond-mat/0010377 (paper H). 
D. Belitz and T. R. Kirkpatrick, Phys. Rev. B 56, 6513 
(1997); D. Belitz, F. Evers, and T. R. Kirkpatrick, Phys. 
Rev. B 58, 9710 (1998). 

See, e.g., J.W. Negele and H. Orland, Quantum Many- 
Particle Systems, Addison Wesley (New York 1988). 
It is important to keep in mind that microscopic details 
negelected in our model can prevent the ferromagnetic 
transition from occurring altogether. For instance, nesting 
of the Fermi surface can lead to antiferromagnetic behavior 
instead. However, we expect that if a ferromagnetic transi- 
tion occurs in a realistic system, then it will be in the same 
universality class as our simple model. 
For the spin-triplet interaction, this approximation is stan- 
dard, and easy to justify. In the spin-singlet channel, 
screening in a metallic system also ensures a short-ranged 
eflFective interaction. As in Ref. ^, we assume that the dis- 
order is not so strong as to put the system in the vicinity of 
a metal-insulator transition, which would cause screening 
to break down. 

For a pedagogical explanation of the replica trick, see, e.g., 
G. Grinstein in Fundamental Problems in Statistical Me- 
chanics VI, edited by E. G. D. Cohen, North Holland (Am- 
sterdam 1985). 

The method employed in Ref. ^ was in turn adapted from 
L. Schafer and F. Wegner, Z. Phys. 38, 113 (1980). 
The physically interesting limit is of course A'' ^ 0, n ^ 00. 
The properties of the homogeneous space in this limit have 
not been investigated. 
^" F. Wegner, Z. Phys. B 35, 207 (1979). 

More precisely, this is the fermionic version of the matrix 
NLo-M, K.B. Efetov, A.I. Larkin, and D.E. Khmelnitskii, 
Zh. Eksp. Teor. Fiz. 79, 1120 (1980) [Sov. Phys. JETP 52, 
568 (1980)]; A.M. Finkel'stein, Zh. Eksp. Teor. Fiz. 84, 168 
(1983) [Sov. Phys. JETP 57, 97 (1983)], while Wegner's 
original model was a bosonic formulation that is adequate 
for noninteracting electrons. 

C. Castellani and C. Di Castro, Phys. Rev. B 34, 5935 
(1986). 

It is possible to integrate out AP in Gaussian approxima- 
tion in closed form while keeping q to all orders. The result, 
however, is not very elucidating, and we therefore expand 
in powers of q before integrating out the massive fiuctua- 
tions. 



16 



■^^ There is also a term proportional to M AP q. However, its 
frequency structure is such that it is effectively proportional 
to the frequency carried by the M-field. As a result, this 
term is less important for our purposes than the M AP 
term, and we neglect it. 

■^^ It is furthermore expected that G and G4, and H and H4, 
as well as all higher coupling constants in the expansion 
of the nonlinear sigma model in powers of q, renormalize 
in the same way. For the sigma model in the absence of 
interactions this is known to be the case (see, e.g., J. Zinn- 
Justin, Quantum Field Theory and Critical Phenomena, 
Clarendon, Oxford 1989), and for the sigma model with 
point- like interactions it is very likely to be true (see, e.g., 
Ref. Since the present model can be reformulated as an 
interacting sigma model, albeit with a non-point like triplet 
interaction amplitude, by integrating out the field M, the 
same property is expected to hold. 

Here we use Ma's method for identifying simple fixed 
points. Accordingly, we use physical arguments to deter- 
mine which coupling constants should be marginal, and 
then check whether this choice leads self-consistently to a 
stable fixed point. See S.-K. Ma, Modern Theory of Critical 
Phenomena (Benjamin, Reading, MA, 1976). 
2^A.B. Harris, J. Phys. C 7, 1671 (1974); J.L. Chayes, L. 
Chayes, D.S. Fisher, and T. Spencer, Phys. Rev. Lett. 57, 
2999 (1986). 

We note that one also obtains this fixed point within a 
theory where the bare value of ad-2 is zero, if one chooses 
ci and C2 to be marginal. This is an interesting alternative 
since it forces one to introduce two different time scales, 
otherwise ci and C2 cannot both be marginal. The problem 
with the formulation of the theory in Ref. ^ was not so 
much the nonlocality of the theory per se as the fact that 
this nonlocality obscured the existence of two time scales. 
■^^ Strictly speaking, SA^^\ Eq. (2.18), is part of the effective 
action. However, none of the quantities we are interested 
in are given in terms of correlations of AP or AA, and the 
couplings between these fluctuations and M or q are all 
irrelevant. We therefore neglect this part of the effective 
action. 

^° R. Shankar, Rev. Mod. Phys. 66, 129 (1994). 
'^^ The arguments given in Ref. ^ for an absence of long-range 
magnetic order were incorrect. 

See, e.g., M.E. Fisher, in Advanced Course on Critical Phe- 
nomena, edited by F.W. Hahne (Springer, New York 1983). 
T.R. Kirkpatrick and D. Belitz, J. Stat. Phys. 87, 1307 
(1997). 

F. Wegner, in Phase Transitions and Critical Phenomena, 
vol.6, edited by C. Domb and M.S. Green (Academic, New 
York 1976). 



17 



